JOHANNES KEPLER UNIVERSITY LINZ Institute of Computational Mathematics A Multigrid Fourier Analysis of a Multigrid Method using Symbolic Computation
نویسندگان
چکیده
For iterative solvers, besides convergence proofs that may state qualitative results for some classes of problems, straight-forward methods to compute (bounds for) convergence rates are of particular interest. A widelyused straight-forward method to analyze the convergence of numerical methods for solving discretized systems of partial differential equations (PDEs) is local Fourier analysis (or local mode analysis). The rates that can be computed with local Fourier analysis are typically the supremum of some rational function. In the past this supremum was merely approximated numerically by interpolation. We are interested in resolving the supremum exactly using a standard tool from symbolic computation: cylindrical algebraic decomposition (CAD). In this paper we work out the details of this symbolic local Fourier analysis for a multigrid solver applied to a PDE-constrained optimization problem.
منابع مشابه
JOHANNES KEPLER UNIVERSITY LINZ Institute of Computational Mathematics Analysis of a Time Multigrid Algorithm for DG-Discretizations in Time
We present and analyze for a scalar linear evolution model problem a time multigrid algorithm for DG-discretizations in time. We derive asymptotically optimized parameters for the smoother, and also an asymptotically sharp convergence estimate for the two grid cycle. Our results hold for any A-stable time stepping scheme and represent the core component for space-time multigrid methods for para...
متن کاملJOHANNES KEPLER UNIVERSITY LINZ Institute of Computational Mathematics A Robust Multigrid Method for Isogeometric Analysis using Boundary Correction
The fast solution of linear systems arising from an isogeometric discretization of a partial differential equation is of great importance for the practical use of Isogeometric Analysis. For classical finite element discretizations, multigrid methods are well known to be fast solvers showing optimal convergence behavior. However, if a geometric multigrid solver is naively applied to a linear sys...
متن کاملEigenvalue Problems in Surface Acoustic Wave Filter Simulations
1 FWF-Start Project Y-192 “3D hp-Finite Elements”, Johannes Kepler University Linz, Altenbergerstraße 69, 4040 Linz, Austria [email protected] 2 Radon Institute for Computational and Applied Mathematics (RICAM), Altenbergerstraße 69, 4040 Linz, Austria [email protected] 3 Institute of Computational Mathematics, Johannes Kepler University Linz, Altenbergerstraße 69, 4040 Linz, Austria ulanger...
متن کاملJOHANNES KEPLER UNIVERSITY LINZ Institute of Computational Mathematics Weighted Poincaré inequalities
Poincaré type inequalities are a key tool in the analysis of partial differential equations. They play a particularly central role in the analysis of domain decomposition and multilevel iterative methods for second-order elliptic problems. When the diffusion coefficient varies within a subdomain or within a coarse grid element, then condition number bounds for these methods based on standard Po...
متن کاملSymbolic Analysis for Boundary Problems: From Rewriting to Parametrized Gröbner Bases
We review our algebraic framework for linear boundary problems (concentrating on ordinary differential equations). Its starting point is an appropriate algebraization of the domain of functions, which we have named integro-differential algebras. The algebraic treatment of boundary problems brings up two new algebraic structures whose symbolic representation and computational realization is base...
متن کامل